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(54) Title: ON-LINE CALIBRATION PROCESS 

^ (57) Abstract: Method for automatic on-line calibration of process models for real-time prediction ol process quality lrom raw 
process measurements by collecting raw process data, processing data collected through a mathematical model to obtain a predic- 
I/") tion of the quality, processing this prediction through two independent dynamic transfers functions thus creating two intermediate 
O signals, storing the two intermediate signals obtained as a function of lime in history, retrieving, at the lime of a real and validated 
^ measurement of the quality, from the history the absolute minimum and maximum values of the two intermediate signals in the lime 
^£ period corresponding to a minimum and maximum specified deadtime, which values define the minimum and maximum prediction 
possible, calculating the deviation as being the difference between the real and validated measurement and the area encompassed 
between the minimum and maximum prediction possible as obtained, and repealing these steps if the absolute value of the deviation 
obtained is zero, or, if the absolute value of the deviation obtained is larger than zero, incorporating the deviation into the process 
^ model and repeating the steps. By using a Kalman filler method for incorporating the deviation into the mathematical model its linear 
^ parameters will be updated, thereby improving the model. The calibration process with the Kalman filler can be applied under non 
steady-state conditions. 
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ON-LINE CALIBRATION PROCESS 



The present invention relates to an automatic on-line 
calibration of input/output models. 

. It is known in process control systems to use one or 
more so-called Quality Estimators (herein after referred 
to as QE) in the real time prediction of certain, 
preferably, key process quality and/or property para- 
meters (normally referred to as outputs) from readily 
available raw process measurements (normally referred to 
as inputs) . Thus, QE is in essence a mathematical 
input/output process designed at predicting relevant 
process values. 

QE ! s are usually identified from collected process 
data. In order to have a useful meaning in real time 
implementation a QE has to be calibrated using historic 
quality measurements, which can be taken on-line or off- 
line depending on the type of process and/or the type of 
measurement envisaged, so as to minimise, or preferably 
avoid, any drift in the predicted quality. QE are 
preferably used in situations, which allow rather 
infrequent and/or delayed measurements of product 
quality. This may be the case when, for instance, the 
amount of time needed to produce the measured value is 
rather long or when the method is relatively costly. 

There are a number of difficulties to be faced in the 
process of automatic on-line calibration of QE such as 
the occurrence of varying or uncertain process/ 
measurement deadtimes and dynamics between the QE inputs 
and the measured qualities as well as a phenomenon 
normally referred to as changing of the process gains, 
i.e. a drift in the ratio between inputs and outputs. 
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In order to combat these unwanted situations, it is 
customary to calibrate QE when the process for which they 
are applicable is in its so-called steady-state, i.e. in 
the situation in which the process fluid is uniform and 
5 constant in composition, state and velocity at the 

entrance and at the exit of the operation. Although such 
calibration will give good results with respect to the 
system to be monitored, it is still considered to be sub- 
optimal as dynamic information available is not used. 

10 This because calibration has to wait until the process 

has reached a steady operating point. Moreover the 
presence of a steady-state detector is required in order 
to know when calibration can start. 

It has now been found that the disadvantages referred 

15 to herein before can be minimized or even overcome by 

applying the process according to the present invention 
which allows for a real time method for automatic on-line 
calibration in a robust manner. The Robust Quality 
Estimator (RQE) according to the present invention 

20 provides a more accurate and robust quality prediction, 

which improves the performance of any quality control 
scheme in which it is applied. For instance, it improves 
the performance of a linear model predictive controller 
when the process is such that the steady-state gains 

25 and/or the dynamics (such as the deadtime) between the 

manipulated variables and the controlled quality are 
varying in an unpredictable manner within certain 
identified boundaries. Moreover, it can also be used to 
facilitate closed-loop control of any process variable 

30 with a difficult dynamic behaviour. 

The present invention therefore relates to a method 
for automatic on-line calibration of a process model for 
real-time prediction of process quality from raw process 
measurements which method comprises: 

35 a) collecting raw process data, 
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b) processing data collected in step a) through the 
process model to obtain a prediction of the quality, 

c) processing this prediction through dynamic transfer 
functions thus creating two intermediate signals, 

5 d) storing the two intermediate signals obtained in 

step c) as a function of time in history, 

e) retrieving at the time of a real and validated 
measurement of the quality from the history the absolute 
minimum and maximum value of the two intermediate signals 

10 in the time period corresponding to a minimum and maximum 

specified deadtime, which values define the minimum and 
maximum prediction possible , 

f) calculating the deviation as being the difference 
between the real and validated measurement and the area 

15 encompassed between the minimum and maximum prediction 

possible as obtained in step e) , 

g) proceeding with step i) if the absolute value of the 
deviation obtained in step f) is zero, or, proceeding 
with step h) if the absolute value of the deviation 

20 obtained in step f) is larger than zero, 

h) incorporating the deviation into the process model, 
and 

i) repeating steps a)-h). 

The process model, which is calibrated with the 
25 method of the present invention, is suitably a so-called 

input-output parametric model, which has been obtained 
off-line from history process data and quality 
measurement. Examples of such models are Multiple Linear 
Regression as described in for example, Introduction to 
30 linear regression analysis by Montgomery and Peck, John 

Wiley & Sons, 1992, Linear Dynamic Model (in the Laplace 
transform Domain) as for example described in Linear 
Systems by Keilath, Prentice-Hall, Information & System 
sciences series, 1980 and Radial Bias Function Neural 
35 Network (optionally with Gaussian function) as for 
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example described in T. Poggio and F. Girosi. Network for 
approximation and learning. The Proceedings of the IEEE, 
78 (9) :1485-1497, September 1990. Depending on the nature 
of the process model applied and the type of raw material 
5 data received, those skilled in the art will select the 

type of process model for the quality estimation best 
fitting the perceived goal. 

Use shall be made of the following Figures to 
illustrate the method of the present invention more 
10 clearly. 

Figure 1 is a functional block-scheme illustrating a 
preferred embodiment of the method of the present 
invention. 

Figure 2 is a graph illustrating different terms of 
15 the method of the present invention. 

Figure 3 is a simplified flow chart of a distillation 
column being controlled by means of a quality estimator. 

Figure 4 shows the dynamic response of the minimum 
and maximum value of the expected quality after a step 
20 change of the set point of the process is performed. 

Figure 5 illustrates the improved distillation 
operation before and after an on-line quality estimator 
is incorporated in the control loop. 

Figure 1 shows a process model (1) having input from 
25 raw process data (2) . The process model (1) provides an 

estimated quality (11), which is used as input for 
controller (12), which may control for example a valve 
(not shown) . Figure 1 also shows a module (3) wherein 
steps (c) and (d) are performed. Further shown is a 
30 validation module (5), which validates the real quality 

measurement (4) to obtain a real and validated quality 
measurement (6) . Based on the input from module (3) and 
the real and validated quality measurement (6) a 
deviation is calculated in (7). If the deviation is 
35 greater than zero as described in step (g) the deviation 
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(8) is used for calibration of the process model (1), 
preferably by making use of the Kalman Filter (9) . 

The collection of raw process data (2) in step (a) to 
be used in the method according to the present invention 
5 can be carried out by methods known in the art. It is 

customary in process control technology to measure 
data (2) at a number of points over a period of time. For 
instance, in refining operations, operating parameters 
such as temperature, pressure and flow are normally 

10 measured at frequent intervals, or even in a continuous 

manner and they can be stored and processed in many ways 
as is known to those skilled in the art. 

In order to get a prediction of the quality (11) out 
of the raw process data (2) collected the above referred 

15 to process model (1) will be used in step (b) . Step (b) 

is thus the quality prediction step. 

Step (c) is an essential step in the method for 
automatic on-line calibration. This and further steps 
will also be illustrated by making use of Figure 2. In 

20 these steps the calculation of the minimum, and maximum 

prediction possible at the time of the real and validated 
measurement (s) of the quality is performed. Step (c) is 
suitably performed by applying two dynamic transfer 
functions (so-called uncertain dynamics) to the 

25 prediction of the quality (11) (the undelayed real time), 

thus creating two intermediate signals. Suitably two or 
more independent dynamic transfer functions are applied. 
Dynamic transfer functions are well known tools to one 
skilled in the art and are for example described in 

30 Linear Systems by Keilath, Prentice-Hall, Information & 

System sciences series, 1980. In step (d) these 
intermediate signals (20, 21) are stored as a function of 
time in history. This will result in essence in an 
(uncertainty) area (22) in which the actual process 

35 response should be placed and which will become very 
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narrow when reaching the steady-state situation (23, 24) . 
It is also possible that the uncertainty area (22), in a 
non-steady state situation, is reduced to a line 
corresponding to the event in which the independent 
5 dynamic transfer functions are identical (this situation 

is not shown in Figure 2) . The so-called minimum and 
maximum prediction possible are obtained by calculating 
from the history the absolute minimum (27) and maximum 
values (28) of these two intermediate signals (20, 21) in 

10 the time period corresponding to a minimum (25) and 

maximum (26) specified deadtime. The deadtime is a 
function of the' virtual location of the quality estimator 
relative to the location where the real quality is 
measured, time for the real quality to be measured and 

15 other process conditions, for example flow rate and 

liquid hold-up. The deadtime can be easily determined by 
one skilled in the art. From this input a maximum (26) 
and minimum (25) deadtime is defined representing the 
time period of the process history in which in step (f) 

20 the real and validated quality measurement (29 -> 29') is 

compared with the predicted quality area (22) and the 
specific minimum (27) and maximum (28) possible quality 
values . 

Before reaching the steady-sate situation, the area 
25 (22) can be very wide. The state of the art systems will 

either only calibrate during steady-state or have the 
risk of making a false calibration in case the real and 
validated measurement (s) of the quality is within the 
above mentioned area. The method according to the present 
30 invention, however, is specifically designed to calibrate 

only when the real and validated measurement (s) (29) of 
the quality are outside the uncertainty area (22), thus 
preventing instabilities in closed-loop. Advantageously 
the calibration method according to the present invention 
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can be performed under steady and non-steady state 
conditions . 

In step (e) in the method according to the present 
invention part of the calibration process is carried out 
5 by calculating the deviation (30) (the so-called 

prediction error) as being the distance between the real 
and validated measurement (29') and the area (22) 
encompassed between the minimum (27) and maximum (28) 
prediction possible as obtained from the earlier 
10 calculation. 

The real and later validated measurement (29) of the 
quality can be an on-line or off-line measurement of the 
quality. The quality to be measured can be properties of 
process streams or products, like for example viscosity, 
15 viscosity index, vapour pressure, boiling points, free 

and cloud points, density, octane number and cetane 
number, and compositional information like for example 
sulphur content, aromatics content, benzene content and 
olefin content can also be measured at frequent intervals 
20 or even continuous manner by making use of off-line or 

on-line analytical apparatuses. Such apparatuses may 
measure the property directly making use of for example 
on-line viscometric analysers, GLC or NMR analysers. The 
quality may also be measured indirectly by making use of 
near infrared prediction methods as for example described 
in EP-A-285251 for measuring octane number. 

In step (g) the usefulness for calibration purposes 
of the real and validated measurement of the quality is 
determined. Only measurements (29') of the quality, which 
are outside the uncertainty area (22), can be used for 
calibration of the model. In other words, if the 
calculation of the deviation (30) as described herein 
above shows that the absolute value of the deviation 
obtained is zero, meaning that the validated and real 
measurement of the quality is within the uncertainty area 
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(22) or more precise, between the minimum (27) and 
maximum (28) possible quality values, the deviation (30) 
found will not be used as further input in the 
calibration process but the system will continue by 
5 repeating the steps carried out up till now as there is 

no need to refine the system* If, however, the deviation 
(30) as calculated shows that the absolute value of the 
deviation (30) is larger than zero, as shown in Figure 2, 
the deviation (30) obtained will be incorporated into the 

10 process model in step (h) and the previous steps will be 

repeated (step i) . The net result will be the generation 
of a modified, more precise, predictive process model, 
which will then serve as the basis for further 
modifications depending on the level of deviations being 

15 observed during the course of the calibrating process. 

Preferably step (h) is performed, such that 
incorporation of the deviation (8) into the process model 
(1) is performed with the use of a Kalman filter (9) (See 
Figure 1) . The result of performing step (h) in such a 

20 manner will be that the deviation can be incorporated 

into the process model by adjusting its linear parameters 
thereby updating the prediction band and improving the 
process model. The use of a Kalman filter is well known 
in the art of process control operations. Reference is 

25 made in this respect to "Stochastic Processes and 

Filtering Theory" by Jazwinski (Academic Press, 
Mathematics and Science and Engineering, Vol. 64, 1970). 
Since Kalman filters are in essence optimal stochastic 
filters they also filter out, or even eliminate, the 

30 noise on the measured quality, which makes them very 

suitable for use in the method according to the present 
invention. 

It should be noted that the use of Kalman filters is 
not limited to calibration operations, which are carried 
35 out under non steady-state conditions, as it is equally 
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capable of providing useful information when a process is 
being operated under steady-state conditions. 

It has been found that by combining the Kalman filter 
with the process according to the present invention an 
5 even more robust control method is obtained. The use of 

Kalman filter has the additional advantage that it will 
keep improving the accuracy of the quality estimation 
process.. In the event that no real and validated 
measurement of the quality is received, calibration as 

10 defined in steps e, f and g is not carried out. The 

system will repeat steps a-d until a further real and 
validated measurement of the quality is received. 

The calibration process as described in the present 
invention can be extrapolated for robust multivariable 

15 predictive controllers to cover uncertain dynamics in the 

control model for all the transfer functions between the 
manipulated variables and the controlled variables. 

A further advantage of the present method is that as 
a by-product of the calibration process (the steps 

20 described above) performed by the method of the present 

invention, the real time accuracy of the prediction is 
exactly known at any time (via the calculated prediction 
error or deviation (30)). This avoids costly and 
extensive quality measurement campaign for validation 

25 purpose as required with conventional Quality Estimators. 

Examples of applications for the present invention 

are: 

Distillation processes, wherein the quality to be 
estimated by the quality estimator is for example the 

30 quality, for example composition, viscosity, density or 

boiling point, of the fractions obtained in the 
distillation. The raw process data as collected in 
step (a) may be the flow rates, feed and product 
temperature, tray temperatures, overhead temperatures, 

35 system temperatures, reflux ratio, circulation reflux 
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duty and the reboiler duty. Exemplary distillation 
processes are crude distillers, vacuum distillers, light 
ends splitters as for example de-ethanisers, de- 
propanisers, de-butanisers, ethane-ethane splitters, 
5 propene-propane splitters, butene-butane splitters. 

Conversion processes, wherein the quality to be 
estimated by the quality estimator is for example the 
quality, for example composition or yield, viscosity, 
density, boiling point, octane number, cetane number, 

10 melt index of the product as obtained in the conversion. 

The raw process data as collected in step (a) may be the 
flow rates, feed composition, density, catalyst age, or 
temperatures. Exemplary conversion processes are 
polymerisation processes, hydrocracking, fluid catalytic 

15 cracking, hydrotreating, hydrogenation, 

hydrodesulphurisation, (hydro) dewaxing, 
hydroisomerisation, reforming or delayed coking. 

Blending processes, wherein the quality to be 
estimated by the quality estimator may be the octane 

20 number or cetane number, viscosity, viscosity index, 

boiling point or composition. The raw process data as 
collected in step (a) may be temperatures of the 
different feed streams or compositional data of the 
buffer tanks upstream. 

25 The application and advantages of the Robust Quality 

Estimator according to the present invention is 
illustrated by the following more detailed description of 
a control of a simple two-cut splitter as shown in 
Figure 3. This Figure shows a schematic representation of 

30 a distillation column (33) provided with a feed (34) and 

top product outlet conduit (31) and bottom outlet conduit 
(35) . The unit is further equipped with a condenser (36) , 
a overhead drum (38) and a reboiler (37) . Various 
controls are present: a level control (39), a top 

35 temperature control (40), a level control (41) and a 
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reboil control (42) . The real quality of the top product 
(31) is measured by a continuous analyzer (32). This 
quality can be for instance a final boiling point or 
impurity content (amount of heavy product in the top 
5 product) . 

Such a distillation process will act as follows on 
for example a step test on the top temperature by 
changing the set point of the top temperature controller 
(40) acting on the reflux valve (42) . The quality 

10 measured by the analyzer (32) will exhibit some first 

order dynamics with deadtime. The deadtime originates 
from transportation from the top of column (33) to the 
analyzer sampling point (32) . This deadtime depends on 
how much product is made and will vary accordingly 

15 (resulting in an uncertain deadtime phenomenon) . If now 

we look at the dynamic response, slow dynamics will occur 
in case the liquid hold-up of the overhead drum (38) is 
small (corresponding to a low level) and vice versa 
(resulting in an uncertain dynamics phenomenon) . The 

20 dynamic response including the deadtime will thus be 

influenced by the level controller LC (39) on the 
overhead drum (38) . From the above the minimum and 
maximum deadtime can be determined. 

Figure 4 shows possible responses to this test step 

25 on the top temperature controller (40) . On the y-axis the 

quality as measured by analyzer (32) is projected. On the 
x-axis time is projected. On time (tl) the temperature 
step is performed as shown by line (51) . Between a 
minimum (52) deadtime and a maximum (53) deadtime the 

30 quality as measured at (32) will start to change from 

steady state quality (56) to a next steady state quality 
(57) . Between steady state quality measurement (56) and 
(57) the quality as measured at (32) can vary between the 
minimum (55) and maximum (56) quality prediction 

35 possible. In a situation as shown in Figure 4 and 
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irrespective of the type of controller type (e.g. a model 
predictive controller) , it will be very tough for the 
controller to achieve tight quality control when the 
controlled variable is the measured quality from the 
5 analyzer (32) . This because of the varying deadtime and 

uncertain dynamics between the top temperature and the 
quality measured by the analyzer as explained above and 
shown in Figure 4 . 

The above disadvantage is overcome by estimating the 

10 quality at the top section of the column with the Robust 

Quality Estimator (RQE) (43 in Figure 3) which is 
calibrated on-line according to the method of the present 
invention. The predictive process model of this RQE is 
using raw process measurements such as key temperatures 

15 in the column, pressures and flow rates. By controlling 

the process based on the quality estimated by the virtual 
analyzer (RQE) (43) instead of the real analyzer all 
tough phenomena, such as varying deadtime/dynamics due to 
the location of the sampling point and/or phenomena such 

20 as varying gain due to change in operating point, have 

been removed from the control loop. The real-analyzer 
(32) in its turn provides the real and later validated 
measurement for robust calibration of the predictive 
process model according to the present invention. Thus 

25 the RQE (43) will provide the controller with an early 

and accurate prediction of the product quality. 

The above method has been applied in practice to a 
commercial benzene-toluene splitter, in this distillation 
unit benzene and toluene are recovered from an extract by 

30 means of distillation. The production objective is to 

maintain the amount of toluene impurity in the benzene 
below some upper limit. When this upper limit is 
exceeded, the product must be re-processed leading tot a 
reduction in production capacity. The commercial toluene/ 

35 benzene splitter was known to be a difficult column to 
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control. The column had a non-linear response. An on-line 
analyzer such as showed in Figure 3 measured the top 
product quality. Typically, the continuous analyzer took 
several minutes to detect changes in the quality of the 
5 product. The response time of this analyzer was also 

dependent upon the production rate through the column. 
This meant that the analyzer was slower to reflect 
changes in product quality during lower production rates. 
The original Multivariable Predictive Controller <MPC) 

10 was not robust against large changes in dead time and 

process non-linearities . 

The MPC controller was subsequently modified to 
control the predictive quality, i.e. the virtual analyzer 
or analyzer predictor (43) instead of the real analyzer 

15 (32) (see Figure 3) . The RQE allows for an early 

prediction of the Benzene quality without having to wait 
for the transportation delay in the overhead line as well 
as the measurement delay due to the analyzer itself. The 
RQE was periodically calibrated according to the present 

20 invention from the on-line analyzer measurement to adjust 

to any process non-linearity/disturbance. Figure 5 
indicates the improvement in performance - tighter 
control - after RQE was incorporated in the control loop 
at day "155". On the y-axis the zero value refers to the 

25 desired or set-point value. The variation on the y-axis 

is the normalized variation of the analyzer (32) output. 
Figure 5 clearly shows that the variation when the 
f control was performed on the quality (32) alone was much 
greater than after day "155", when the RQE (43) was used 

30 to control the quality of product (31) . In hard figures 

the Aromatics plant had been producing approximately 
640 t/d (200 ktpa) of Benzene. The incorporation of the 
RQE in the control loop has realized a situation where 
the plant now has a capability to produce at least 

35 850 t/d (287 ktpa) . This has led to a record production 
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year, despite a 37-day plant shutdown for major planned 
maintenance work and later in the year having to reduce 
production somewhat due to the lack of available feed 
from outside of the refinery. 
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C L A IMS 



1. A method for automatic on-line calibration of a 
process model for real-time prediction of process quality 
from raw process measurements which method comprises 

a) collecting raw process data, 

b) processing data collected in step a) through the 
process model to obtain a prediction of the quality, 

c) processing this prediction through dynamic transfer 
functions thus creating two intermediate signals, 

d) storing the two intermediate signals obtained in 
step c) as a function of time in history, 

e) retrieving at the time of a real and validated 
measurement of the quality from the history the absolute 
minimum and maximum value of the two intermediate signals 
in the time period corresponding to a minimum and maximum 
specified deadtime, which values define the minimum and 
maximum prediction possible, 

f ) calculating the deviation as being the difference 
between the real and validated measurement and the area 
encompassed between the minimum and maximum prediction 
possible as obtained in step e) , 

g) proceeding with step i) if the absolute value of the 
deviation obtained in step f) is zero, or, proceeding 
with step h) if the absolute value of the deviation 
obtained in step f) is larger than zero, 

h) incorporating the deviation into the process model, 
and 

i) repeating steps a)-h). 

2. A method according to claim 1, in which as process 
model a Multiple Linear Regression model is used. 

3. A method according to claim 1, in which as process 
model a Linear Dynamic Model is used. 
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4. A method according to claim 1, in which as process 
model a Radial Basis Function Neural Network is used. 

5. A method according to one or more of claims 1-4, in 
which in step h) a Kalman filter method is used to 
incorporate the deviation into the process model by 
adjusting its linear parameters thereby upgrading the 
prediction and improving the process model by self 
learning . 

6. A method according to claim 5, in which the Kalman 
filter is used in step h) under non steady-state 
conditions of the process. 

7. A method according to any one of claims 1-6, wherein 
the prediction in step (c) is processed through two 
independent dynamic transfer functions. 
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